The Chemical Composition of Cernis 52 (BD+31 640) 
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We present an abundance analysis of the star Cernis 52 in whose spectrum 
we recently reported the napthalene cation in absorption at 6707.4 A. This star 
is on a line of sight to the Perseus molecular complex. The analysis of high- 
resolution spectra using a ^-minimization procedure and a grid of synthetic 
spectra provides the stellar parameters and the abundances of O, Mg, Si, S, Ca, 
and Fe. The stellar parameters of this star are found to be T c g = 8350 ± 200 K, 
log(#/cm s 2 )= 4.2 ± 0.4 dex. We derived a metallicity of [Fe/H] = -0.01 ± 0.15. 
These stellar parameters are consistent with a star of ~ 2 M Q in a pre-main- 
sequence evolutionary stage. The stellar spectrum is significantly veiled in the 
spectral range AA5150 — 6730 A up to almost 55 per cent of the total flux at 
5150 A and decreasing towards longer wavelengths. Using Johnson-Cousins and 
2MASS photometric data, we determine a distance to Cernis 52 of 231+gg 5 pc 
considering the error bars of the stellar parameters. This determination places 
the star at a similar distance to the young cluster IC 348. This together with 
its radial velocity, v r = 13.7 ± 1 kms -1 , its proper motion and probable young 
age support Cernis 52 as a likely member of IC 348. We determine a rotational 
velocity of v sin % = 65 ± 5 km s _1 for this star. We confirm that the stellar 
resonance line of Li I at 6707.8 A is unable to fit the broad feature at 6707.4 A. 
This feature should have a interstellar origin and could possibly form in the dark 
cloud L1470 surrounding all the cluster IC 348 at about the same distance. 

Subject headings: stars: abundances — individual: [C93]52 — individual: BD+31 640 
— ISM:molecules — ISM:lines and bands — ISM:abundances 
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1. Introduction 

In a recent study of the diffuse inter- 
stellar bands (DIBs) toward a region of 
anomalous microwave emission, Iglesias- 
Groth et al. (2008) have shown evidence 
for the presence of the naphthalene cation, 
Ci Hg , the simplest Policyclic Aromatic 
Hydrocarbon (PAH), in the line of sight 
towards the moderately reddened star Cer- 
nis 52 (BD+31° 640, V = 11.4, E(B-V) = 
0.9, Cernis 1993). The PAHs, including 
the naphthalene cation, have been also de- 
tected in a cometary dust sample returned 
to Earth by Stardust and in interplanetary 
dust particles, possibly of cometary origin 
(see the review by Li 2009). 

Cernis 52 is an early type star, classi- 
fied as A3V by Cernis (1993), has coordi- 
nates a = 03 h 43 m 00.3 s and 5 = +31 58'26" 
(J2000.0), and is located at an angular 
separation of less than one degree from 
the very young stellar cluster IC 348 in 
the Perseus OB2 molecular complex. The 
photometric distance to star Cernis 52 is 
240 pc (Cernis 1993), consistent with a 
location in the molecular complex OB2 
where the microwave emitting cloud is 
also most likely located (see Watson et 
al. 2005). The uncertainties prevent to 
conclude whether the star is embedded in 
a cloud of this young star forming region 
or lay behind. 
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Iglesias-Groth et al. (2008) claimed that 
the relative strength of several diffuse in- 
terstellar bands detected in the spectrum 
of Cernis 52 was consistent with the pres- 
ence of a molecular cloud in the line of 
sight. The anomalous microwave emission 
detected in the direction to Perseus (Wat- 
son et al. 2005, 2006) could be associated 
to electric dipole radiation of fast spinning 
hydrogenated carbon-based molecules in a 
molecular cloud, a mechanism originally 
proposed by Draine and Lazarian (1998). 
It is important to establish whether the 
cloud responsible for the excess color of 
Cernis 52 hosts the carriers of both, the 
diffuse interstellar bands observed in the 
spectrum of this star and the anomalous 
microwave emission detected in its line of 
sight. 

Iglesias-Groth et al. (2008) also noted 
the near-coincidence in wavelength be- 
tween the interstellar cation's feature and 
the stellar Li I resonance line. To gain 
insight into the contribution of the stellar 
line to the observed broad feature, we un- 
dertook a thorough analysis by spectrum 
synthesis of Cernis 52's spectrum with two 
principal goals - a general abundance anal- 
ysis for the star and spectrum synthesis 
fits to the spectrum around 6707 A in or- 
der assess the Li I line's contribution. 

Here we present a precise determination 
of the stellar parameters and a detailed 
chemical composition study of the star 
Cernis 52. This abundance analysis should 
provide indications of the likely Li abun- 
dance for the star, i.e., if the star is not 
chemically peculiar should have a lithium 
abundance no greater than the cosmic Li 
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abundance of A(Li]H= 3.3. In this paper 
we will discuss with special attention the 
naphthalene cation's feature at 6707.4 A 
in the spectrum of Cernis 52 and investi- 
gate the relationship of this star with the 
Perseus star forming region. 

2. Observations 

We have made use of spectroscopic data 
obtained with the 2dcoude cross-dispersed 
echelle (CS23, Tull et al. 1995) at the 2.7 
m Harlan J. Smith Telescope at McDonald 
Observatory (Texas, USA). We observed 
the star Cernis 52 on 26, 27 August 2007, 
and 1, 2 March 2008 in the spectral range 
AA3740 - 10160 A, with resolving power 
X/5X ~ 60, 000 and ~ 30, 000, respectively. 
These two datasets were independently re- 
duced and properly combined to produce 
final spectra of SNR ~150 and 250, respec- 
tively, at the position of Ha for a dispersion 
of 0.05 A/pixel. 

We also performed spectroscopic ob- 
servations with the cross-dispersed echelle 
spectrograph SARG (Gratton et al. 2001) 
at the Telescopio Nazionale Galileo (TNG) 
at the Observatorio del Roque de los 
Muchachos (La Palma, Canary Islands, 
Spain). We observed the target star on 
4 November 2006 in the spectral ranges 
AA4600 - 7975 and AA5600 - 10100 A, with 
a spectral resolution of X/5X ~ 57, 000 and 
a final SNR ~ 100 at Ha for a dispersion 
of 0.034 A/pixel. 

Early-type rapidly rotating hot stars 
were also observed for correction of pos- 
sible telluric absorptions in all campaigns. 
All the spectra were reduced within IRAF, 
wavelength calibrated and co-added after 



1 A(Li)= log[iV(Li)/iV(H)] + 12 



correction for telluric lines. 

3. Rotational and radial velocity 

Firstly, each of the three final spectra 
was cross-correlated (within the package 
MOLLY) with a synthetic spectrum prop- 
erly broadened with a v sinz = 80 kms -1 . 
Thus, we corrected these spectra for their 
radial velocities and derived a rotational 
velocity of v sini = 65 ± 5 km s -1 by com- 
paring broadened and veiled versions of 
a A3V synthetic spectrum in steps of 
5 kms -1 , adopting an spherical rotational 
profile with linearized limb-darkening e = 
0.46, appropriate to the spectral type of 
the star (Al-Naimiy 1978). 

Finally, we estimate a radial velocity 
of Cernis 52 at v r = +13.7 ± 1 kms -1 . 
The cross- correlation was performed in the 
spectral range AA5140 - 6700 A and mask- 
ing the broad Ha profile and the inter- 
stellar features at 5889-95 A and 6280 A. 
This determination is consistent with the 
mean radial velocity of the young cluster 
IC 348, which has been recently estimated 
at 13.5 ± 1.2 kms" 1 from 30 near solar- 
mass cluster members (Dahm 2008). How- 
ever, Nordhagen et al. (2006) derived a 
mean cluster velocity of 15.9 ± 0.8 km s -1 , 
also from 30 cluster members. 

We also determine a radial velocity of 
14.2 ± 1.4 kms -1 from the narrow atomic 
interstellar lines Na I D 5889-95 A, Li I 
6708 A and K I 7698 A. The uncertainty 
on the velocity comes from the dispersion 
of these velocity measurements. 

These determinations share the veloc- 
ity of the most luminous cluster member 
HD281159 (~ 14 kms" 1 , Evans 1967) and 
an expanding supernova shell in which the 
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young cluster IC 348 is embedded (~ 12 — 
15 kms -1 , Snow et al. 1994). In addition, 
these atomic interstellar lines have a veloc- 
ity consistent with that of the interstellar 
hyperfine transitions of the molecule NH 3 
measured in radio wavelegths (Rosolowsky 
et al. 2008): the LSR velocities are 7.5 ± 
1.4 kms -1 from the ISM lines and 6.1 — 
8.4 km s _1 from the radio transitions. 




4. Companion star 

Using the camera FastCam (Oscoz et 
al. 2008) attached to the Nordic Optical 
Telescope (NOT) at the Observatorio del 
Roque de los Muchachos, we obtained im- 
ages in V, R and / bands of Cernis 52 on 
25 October 2008. We found a companion 
star at a distance of 818 ± 7 mas. This 
camera is able to obtain and record series 
of thousands of short exposures, and with a 
specially designed software can reach angu- 
lar resolutions of roughly 0.1". In Fig. [I] we 
display the resulting image in the R band. 
We have zoomed in the image to show 
the two stars clearly resolved. The cam- 
era has 512x512 pixels and a pixel size of 
31.01±0.01 mas. The companion star is lo- 
cated in the North- West of Cernis 52, with 
a position angle North-East 9 = 325.9°. 

The companion star is roughly 1.7 ±0.2 
mag fainter than Cernis 52 in the R band. 
In our spectroscopic observations with the 
instrument CS23, obtained prior to the 
discovery of this companion star, the slit 
widths were 1.2" and 2.4" for the 2007 
and 2008 campaigns, respectively. The slit 
length was 8.2" in both campaigns. Since 
the slit width was quite large and the see- 
ing during the observations was typically 
~ 2 — 2.5", most of the light from the 
companion star must have entered in the 
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Fig. 1. — FastCam image in the R band 
of Cernis [93] 52. The displayed field size is 
~ 3.5" x 3.5". There exists a companion 
star at a distance of 818 ± 7 mas in the 
direction North- West of Cernis 52. 

slit and contaminated the spectrum of Cer- 
nis 52. Assuming that the primary is A3V 
and that both stars are a physical pair we 
can estimate the spectral type of the com- 
panion from the observed magnitude dif- 
ference in the R-band. We adopt a T c s = 
8500 K and log g = 4.0 for the primary. 
Thus, we find a probable T c s = 5800 K 
for the companion from suitable bolomet- 
ric corrections (Bessell et al. 1998) and 
assuming that both stars share the same 
color excess. Thus, we search for evidence 
of the spectrum of a solar-type compan- 
ion in the 2007 and 2008 final spectra of 
Cernis 52. We cross-correlated these spec- 
tra with the spectrum of the Sun (Kurucz 
et al. 1984), properly broadened with dif- 
ferent rotational velocities and in a radial 
velocity range from -150 to 150 km s -1 , but 
we did not find any clear signature of this 
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star in the cross-correlation function. 

It is likely that the companion star have 
slightly contaminated the spectrum of Cer- 
nis 52. From the magnitude difference 
in the R band and taking into account 
that the companion star was partially out 
from the slit, we estimate the flux ratio of 
both stars in our spectrum at F sec /F C52 = 
0.11 ± 0.05 in the wavelength region corre- 
sponding to the R band. However, as we 
will see this is not sufficient to explain the 
level of veiling found in the analysis of the 
photospheric features of Cernis 52. There 
must be another reason for the apparent 
weakness of the stellar features. 

In Fig. [2] we display an optical image 
from the Palomar Digital Sky SurvejH of 
Cernis 52 with a field size of 5' x 5'. In this 
DSS-2 red image one can see an entended 
emission of interstellar gas surrounding the 
star which may be the reason why the spec- 
trum of Cernis 52 is veiled. Our spectra 
were not background subtracted, since the 
star is very bright and the raw images do 
not show any clear sky line. However, in 
the raw images there is neither any indica- 
tion of the presence of this nebulosity, ap- 
part from the high number of counts of the 
sky, and any signal must have been sub- 
tracted when correcting for the scattered 
light. In Fig. [2] we see that the interstellar 
cloud is more prominent towards the North 
of the star, ~ 0.5' from the star center. 
However, we do not see any asymmetry in 
the signal of the sky spectrum which may 
indicate that the slit was not positioned 
South— North. The interstellar cloud may 
have contaminated our spectra but we sus- 
pect that most of the light coming from the 

2 http: / /archive. eso.org/ dss / dss. 



cloud must have been subtracted during 
the data reduction. A discussion of the in- 
terstellar emission in the proximity of Cer- 
nis 52, using Spitzer data, is deferred to a 
forthcoming paper (S. Iglesias-Groth et al. 
2009, in preparation). 




Fig. 2. — DDS-2 image in the R band of 
Cernis [93]52. The displayed field size is 
~ 5' x 5'. 

In the next section, we will concentrate 
on the determination of the stellar param- 
eters of Cernis 52, taking into account any 
possible source of veiling in the spectrum, 
whose flux we will denote as F vea . 

5. Stellar Parameters 

We have inspected the spectrum of Cer- 
nis 52 to search for atomic stellar lines 
appropriate to provide accurate element 
abundances. We realized that the stel- 
lar lines were significantly veiled even at 
a level of ~ 50 per cent at 5180 A where 
the Mg lb triplet is located. Thus, the Fe 
abundances obtained from the lines Fe I 



5 



5214 A and Fe II 5273 A if we do not in- 
clude any veiling are [Fe/H] = —0.55 and - 
0.8, respectively. This metallicity seems to 
be too low for the location of this star in a 
star-forming region. In addition, the Mg lb 
5172 A would require an abundance of 
[Mg/Fe]~ —1.0. These results are achieved 
for all reasonable values of the atmospheric 
parameters, T e s and log g. The analy- 
sis of this spectrum requires a tool able to 
take into account the effect of any possible 
veiling when computing atmospheric abun- 
dances. To simplify, we assumed the veil- 
ing to be a linear function of wavelength, 
and thus defined by two parameters, its 
value at 6562 A, / 6562 = F^ 2 /F^f, and 
the slope, tuq. Note that the total flux is 
defined as F to tai = ^Veii + ^C52, where 

-^veil 

and Fc52 are the flux contributions of the 
source of veiling and the continuum of Cer- 
nis 52, respectively. We decided to split 
the determination of the stellar and veil- 
ing parameters in two steps: 1) using the 
Ha profile to derive the stellar tempera- 
ture and the veiling at 6562 A; and 2) us- 
ing Fe I and Fe II features to derive the 
surface gravity, the slope of the veiling and 
the metallicity of Cernis 52. 



Table 1: Ranges and steps of model param- 
eters 



Parameter 


Range 


Step 


T eff (K) 


7850 - 


-> 8750 


100 


log(g/cm s 2 ) 


3. - 


-> 5 


0.1 


[Fe/H] 


-0.5 


-> 0.5 


0.05 


^6562 


-» 


1.0 


0.05 


m (A- 1 ) 


-> - 


-0.001 


-0.0001 



1.10 


'± i . i ii n 


.i 




i. 


1.00 








0.90 








0.80 


T„= 

: IV 


= 8050K %i 
-- 8350K X 




0.70 


r — I 

H i , , , 


= 8650K 1 j 


Cernis 52 J 



6540 6550 6560 6570 6580 6590 
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Fig. 3. — Synthetic Ha profiles for three 
effective temperatures compared to the ob- 
served spectrum of Cernis 52 (SNR ~ 150) 
and normalized to the level of the observed 
spectrum of Cernis 52 at 6540 A. We also 
depict the fitting regions at the top. 

5.1. Ha Profile 

The wings of Ha are a very good tem- 
perature indicator (e.g., Barklem et al. 
2002). Adopting the theory of Ali & Griem 
(1965, 1966) for resonance broadening and 
Griem (1960) for Stark broadening, we 
computed Ha profiles for several effective 
temperatures, using the code SYNTHE 
(Kurucz 2005; Sbordone 2005). For further 
details on the computations of hydrogen 
lines in SYNTHE, see Castelli & Kurucz 
(2001) and Cowley & Castelli (2002). 

In this section, we will use the spec- 
trum of Cernis 52 taken in 2007 since the 
2008 spectrum had a different wavelength 
setting and the Ha line appears in an in- 
terorder region. In Fig. [3] we compare the 
synthetic Ha profiles with the observed 
profile for several temperatures. To eval- 
uate the goodness of fit, we employ a re- 
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duced %2 statistic, 

i=l x ' 

where A" is the number of wavelength 
points, M is the number of free parameters 
(here two, T e fj and fesw), fi is the synthetic 
normalized flux, Fi is the observed normal- 
ized flux, and <7; = 1/SNR. The SNR was 
estimated as a constant average value in 
continuum regions close to the observed 
Ha profile. The fitting regions are indi- 
cated in Fig. [31 which contains all of the 
spectral regions close the center of the Ha 
profile where there are no stellar lines and 
where the normalized flux is greater than 
~ 0.7. This figure also shows that the sen- 
sitivity of this method to the T e g decreases 
toward higher effective temperatures. 

In order to estimate the uncertainties on 
temperature and veiling, we built 1000 re- 
alizations of the observed spectrum taking 
into account the SNR. We used a boost- 
rap Monte Carlo method to define the 1- 
a confidence regions and found as most 
likely values T eff = 8300 ± 100 K and 
^6562 = 0.50 ± 0.05. The corresponding 
histograms are displayed in Fig. HI The 
T e Q and /6562 were varied within the ranges 
given in Table [TJ To check the influence 
of the spectrum of the companion star we 
perform another simulation by introduc- 
ing the Ha profile of a pre-main sequence 
solar metallicity star with T e fj = 5800 K 
and log g = 3.8, and assuming that this 
spectrum was contributing with the 10% 
of the stellar flux. The adopted surface 
gravity for the companion star coincides 
with the expected value for a pre-main se- 
quence star of this temperature (see Sec- 
tion [6]). We use a rotational velocity of 



vsini = 100 kms -1 and a limb-darkening 
coefficient e = 0.6, because the spectrum 
of Cernis 52 does not show any stellar nar- 
row absorption and also to be consistent 
with the synthetic spectra computed for 
Section [61 However, we check that us- 
ing lower rotational velocities does not af- 
fect the derived T eS . The result is shown 
in Fig. The distribution of temper- 
atures slightly changes but our adopted 
effective temperature remains almost the 
same, T e g = 8350 ± 100 K. In contrast, 
the required veiling factor changes, being 
in this case / 6562 = 0.40 ± 0.05. 

5.2. Ionization Equilibrium of Iron 

Once we have determinated the veil- 
ing at 6562 A and the effective tempera- 
ture we try to derive the surface gravity 
and metallicity using Fe I and Fe II fea- 
tures of the spectrum of Cernis 52. In 
this case, we will neglect the presence of 
a companion spectrum since it hardly af- 
fects the spectrum of Cernis 52. We use 
a technique which compares a grid of syn- 
thetic spectra with the observed spectrum, 
via a ^-minimization procedure. In this 
case, we will use the spectrum of Cernis 52 
taken in 2008, since it has higher SNR. 
This procedure takes into account the ef- 
fect that veiling causes on the stellar lines 
and it has been applied in X-ray binaries 
to derive the stellar parameters of the com- 
panion stars in these systems (A0620-00, 
Gonzalez Hernandez et al. 2004; Centau- 
rus X-4, Gonzalez Hernandez et al. 2005; 
XTE J1118+480, Gonzalez Hernandez et 
al. 2006, 2008b; and Nova Scorpii 1994, 
Gonzalez Hernandez et al. 2008a). How- 
ever, we needed to slightly modify this pro- 
gram in order to allow a fixed veiling at a 
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given wavelength. In the previous version 
of this program, the slope of the veiling 
was defined using a given veiling at 4500 A. 
We typically produce 10 different values for 
the slope, whose maximum value is zero, 
i.e. a constant veiling, and the minimum 
value provides a veiling of / 6562 = 0.50 and 
/5000 = 2 (see Tabled]). 




8000 8250 8500 8750 
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Fig. 4. — Normalized distributions ob- 
tained for the effective temperature, T e g, 
measured from the Ha profile and the veil- 
ing at 6562 A, ^6562, using Monte Carlo 
simulations. The total number of simula- 
tions was 1000. 

We already know the veiling at 6562 A, 
but we need also to find the slope of the 
veiling in order to properly model the 
observed features. We selected 10 spec- 
tral features of iron in the spectral range 
AA5270 - 6400 A, containing in total 55 
lines of Fe I and 25 lines of Fe II with exci- 
tation potentials between 0.5 and 10.5 eV. 
We compute synthetic spectra with the 
local thermodynamic equilibrium (LTE) 
code MOOG (Sneden 1973), adopting the 
atomic line data from the Vienna Atomic 
Line Database (VALD, Piskunov et al. 
1995) and using a grid of LTE model at- 
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Fig. 5. — Same as Fig. HI but including the 
spectrum of the Sun with a contribution of 
10% of the stellar flux. 

mospheres (Kurucz 1993). The oscillator 
strengths of the iron lines as well as other 
element lines used in this work were ad- 
justed until reproducing both the solar 
atlas of Kurucz et al. (1984) with solar 
abundances (Grevesse et al. 1996) and the 
spectrum of Procyon with its derived abun- 
dances (Andrievsky et al. 1995). We gener- 
ated a grid of synthetic spectra for the Fe I 
and Fe II features varying as free parame- 
ters, the star effective temperature (T e fj), 
the surface gravity (logg), the metallicity 
([Fe/H]), and the veiling, charaterized by 
the / 6562 and m . The microturbulence, 
£ = 2kms _1 , was fixed in each atmo- 
spheric model according to the spectral 
type of Cernis 52 (see e.g. Fossati et al. 
2008). We choose an error of 0.5 km s~ x for 
the microturbulent velocity. 

We fixed the veiling parameter at / 6562 = 
0.50 and allowed the effective temperature 
to be in the range 8250 < T eff < 8450. The 
other three free parameters were varied us- 
ing the ranges and steps given in Table [TJ 
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Fig. 6. — Normalized distributions obtained for the effective temperature, T e g, surface grav- 
ity, log g, metallicity, [Fe/H], and veiling slope, m , using Monte Carlo simulations. The 
labels at the top of each bin indicate the number of simulations consistent with the bin value. 
The total number of simulations was 1000. 



We use a bootstrap Monte-Carlo method 
to derive the resulting histograms of these 
parameters shown in Fig. El correspond- 
ing again to 1000 realizations. We found 
strong evidence of high veiling in the spec- 
tral region under analysis (5270-6400 A). 
This veiling increases toward shorter wave- 
lengths, being as high as / 5 i 8 o ~ 1.2 at the 
Mg lb triplet 5167-83A. This translates 
into a veiling of roughly 54% when com- 
puting the ratio of the flux to the 
total flux, F™£. Therefore, at 5180 A the 
star contributes almost with half of the 
total flux. By inspecting the histograms 
displayed in Fig. [HI we decided to choose 
as most likely values: T eS = 8350 ± 200 



K, log(#/cm s 2 ) = 4.2 ± 0.4, [Fe/H] = 
0.20 ± 0.15, / 6562 = 0.50 ± 0.10, and 
m = -0.00055 ± 0.00010. These stellar 
parameters are consistent with the spectral 
type already reported by Cernis (1993). 

The surface gravity is the stellar param- 
eter that is usually more difficult to deter- 
mine when dealing with fast rotating and 
veiled stars (see e.g. Gonzalez Hernandez 
et al. 2005, 2008b). We warn the reader 
that the error bar given for the surface 
gravity comes from the range of values that 
shows the histogram of this parameter in 
Fig El 
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Table 2 

Chemical Abundances of Cernis 52 



Species 


log £ (X) Q a 


[X/H] 


[X/Fe] 


<7 


A a 


A T off 


A log 9 


A vc ;i 


A « 


A [X/H] 


N h 


Ol 


8.74 


0.54 


0.55 


0.28 


0.20 


0.00 


0.03 


0.15 


-0.10 


0.27 


2 


Mgl 


7.58 


0.35 


0.36 


0.25 


0.25 


0.20 


-0.05 


0.15 


-0.30 


0.47 


5173 


Si II 


7.55 


-0.10 


-0.09 


0.28 


0.20 


0.00 


0.03 


0.10 


-0.15 


0.27 


2 


Si 


7.33 


-0.19 


-0.18 


0.04 


0.02 


0.10 


0.00 


0.05 


0.00 


0.11 


2 


Cal 


6.36 


-0.02 


-0.01 


0.22 


0.08 


0.17 


-0.02 


0.07 


-0.05 


0.21 


7 


Ca 11 


6.36 


0.10 


0.11 


0.25 


0.25 


0.10 


0.10 


0.05 


0.00 


0.29 


5307 


Fe I 


7.50 


0.03 


0.04 


0.19 


0.05 


0.13 


0.01 


0.05 


-0.03 


0.15 


13 


Fell 


7.50 


-0.06 


-0.05 


0.08 


0.03 


0.07 


0.04 


0.04 


-0.06 


0.11 
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Note. — Chemical abundances of Cernis 52 and uncertainties produced for At h = +200 K, Ai og9 = +0.4 
dcx, A vc ii = +0.1, and Aj = +0.5 kms" 1 . All the element abundances were determined by fitting the 
observed spectra with synthetic spectra computed with the LTE code MOOG. 

a The solar element abundances were adopted from Grevesse et al. (1996) and Ecuvillon et al. (2006) 

b Number of spectral features of this atomic specie in the star, or if there is only one, its wavelength. 



6. Photospheric abundances 

Using the derived stellar and veiling pa- 
rameters, we firstly determined the abun- 
dance of Fe from each individual feature of 
Fe I and Fe II in the observed spectrum 
(see Table [2]) . We made use of a y2 pro- 
cedure to find the best fit to the observed 
features as in Gonzalez Hernandez et al. 
(2004, 2008a). In Fig. [7] we show some 
of the spectral regions analysed to obtain 
the Fe abundance. To determine the abun- 
dances of other elements, we adopted the 
average abundance of the Fe abundances 
extracted from Fe I and Fe II features, 
[Fe/H] = -0.01 ± 0.15, for the metallicity 
of the model. 

In Table |2] we provide the average abun- 
dance of each element together with the er- 
rors. The errors in the element abundances 
show their sensitivity to the uncertain- 
ties in the effective temperature (At cS ), 
surface gravity (A log3 ), microturbulence 
(A^), veiling (A vcil ) and the dispersion of 



the abundance measurements from differ- 
ent spectral features (A CT ). The errors A CT 
were estimated as A CT = a/y/N, where a is 
the standard deviation of the N measure- 
ments. For those element species for which 
only one spectral feature was available we 
adopt a = 0.25. The total error in the 
abundance of each element is determined 
by adding in quadrature all these individ- 
ual errors. 

In Figs. El M and M we display several 
spectral regions of the observed spectrum 
of Cernis 52 in comparison with synthetic 
spectra computed using the derived abun- 
dances, except for oxygen, for which the 
best fit abundance was used for each fea- 
ture. In these figures, we also show the 
combined synthetic spectrum of Cernis 52 
and the companion star in which the com- 
panion star only contributes with ~ 10% 
of the stellar flux, and using the same veil- 
ing factors as for the single synthetic spec- 
trum of Cernis 52. The stellar parameters, 
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Fig. 7. — Best synthetic spectral fits (solid lines) to the observed spectrum of Cernis 52 for 
two spectral regions. We also display the combined spectrum (dashed line) including the 
spectrum of Cernis 52 contributing with 92% (top panel) and 91% (bottom panel) of the 
stellar flux and the spectrum of the companion star contributing at 8 and 9%, respectively. 
Note that the veiling factors / 5 i 80 = 1.25 and /5370 = 1.15 have been adopted. The label 
"DIB" refers to diffuse interstellar bands, whereas other labels show stellar absorption lines 
of Cernis 52. 



metallicity and radial velocity of the com- 
panion star were already specified in Sec- 
tion 15.11 We assumed that the radial ve- 
locity of the companion star is the same as 
for Cernis 52. For most of the features, the 
presence of the companion do not change 
significantly the line profiles, except for the 
Mg lb 5167-83A lines and the Ca I fea- 
tures. We note that we have not taken into 



account the presence of the companion star 
in the abundance determination, since this 
effect would be negligible and in any case, 
the abundances would remain within the 
error bars given in Table [2J 
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Fig. 8. — Same as Fig. [7J but for other spectral regions. In this case, the best fit abundance 
of each feature has been used, being [O/H] = 0.74 in the upper panel and [O/H] = 0.34 in the 
lower panel. We also display the combined spectrum (dashed line) including the spectrum 
of Cernis 52 contributing with 90% (top panel) and 87% (bottom panel) of the stellar flux 
and the spectrum of the companion star contributing at 10 and 13%, respectively. Note that 
the veiling factors / 6 i 57 = 0.75 and f 7 j 73 = have been adopted. 



7. Discussion 

7.1. Membership in IC 348 

As shown in Section [31 the radial ve- 
locity of Cernis 52 is consistent with the 
mean radial velocity of known members of 
the cluster IC 348. 

The membership of a star in a clus- 
ter is generally established if the star is 



located nearby on the plane of the sky 
and its radial velocity and proper motion 
coincide with those of the cluster. For- 
tunately, the proper motion of Cernis 52 
has been measured with high accuracy by 
Roser et al. (2008), providing (/j, a cosS, 
Us = (+7.22,-8.62) ± (1.5,1.6) mas yr" 1 . 

The proper motion of the cluster IC 348 
has been determined by Scholz et al. 
(1999) in several systems. In particu- 
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Fig. 9. — Same as Fig. [7J but for other spectral regions. We also display the combined 
spectrum (dashed line) including the spectrum of Cernis 52 contributing with 90% (both 
panels) of the stellar flux and the spectrum of the companion star contributing at 10%, 
respectively. In the bottom panel, the solid line show a synthetic spectrum of Cernis 52 with 
A(Li) = 3.3, whereas the dashed line shows a combined spectrum in which both Cernis 52 
and the companion star have A(Li) = 3.3. Note that the veiling factors ^6360 — 0.60 and 
^6720 — 0.40 have been adopted. The label "ISM Lil" gives the location of the atomic Li 
interstellar line. 



lar, using the Hipparcos system, they de- 
rive a mean proper motion of (^cos^, 
H S = (+4.6,-8.3) ± (3.7,1.6) mas yr" 1 , 
from nine Hipparcos stars of the cluster 
IC 348. 

More recently, the measurement of 
proper motion of this cluster has been im- 
proved considerably by Loktin & Beshenov 



(2003). They determine a proper motion of 
(// Q cos5, n s = (+6.87,-9.15) ± (0.56,0.49) 
mas yr -1 . 

The radial velocity and proper motion 
of Cernis 52 strongly support its member- 
ship in the young cluster IC 348. In ad- 
dition, the nine stars studied by Scholz 
et al. (1999) provide a mean distance of 
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Table 3 
Distance to Cernis 52 



Parameter 


£?(B-V) a =0.9 


Teg b = 8550 K 1 


og g b = 3.8 


log g b = 4.6 


E(B-V) b =0.8 


M* b = 2.2 M 


dm v C (PC) 


237 ±9 


247 ±9 


378 ± 14 


149 ± 9 


273 ± 5 


249 ± 10 


dm R (PC) 


229 ± 10 


237 ± 10 


363 ± 16 


144 ± 11 


256 ±6 


240 ±11 


dmj (pc) 


247 ± 10 


253 ± 11 


390 ± 17 


157 ± 11 


268 ±6 


259 ± 11 


dmj (PC) 


233 ±2 


237 ±2 


369 ±3 


147 ±2 


242 ± 1 


245 ± 2 


d m „ (pc) 


232 ±2 


235 ±2 


367 ±4 


147 ±3 


238 ± 1 


243 ± 2 


dm K (pc) 


227 ± 1 


230 ± 1 


359 ±2 


144 ±2 


231 ± 1 


238 ± 1 


4 (PC) 


231 ±8 


238 ±8 


366 ± 12 


146 ± 5 


257 ± 18 


252 ±9 



a The distance is estimated by adopting T c g = 8350 K, log(g/cm s 2 ) = 4.2 dex, M* = 2 M©, and the color excess E(B — V) = 0.9. 
b The distance is estimated by only changing the parameter indicated in the head of each column. 
c The error bar in the distance is estimated from the uncertainty in the magnitude. 

d Avcragc distance weighted by the uncertainty of each individual distance determination for the all filters. The error bar shows the 
standard deviation from the mean of all the individual determinations. 



261^23 P c - Therefore, Cernis 52 must be 
located at this distance. In the following 
subsections we will discuss the mass, ra- 
dius, age, luminosity and distance of Cer- 
nis 52. 

7.2. Mass, Radius and Age 

In Fig. [10] we depict theoretical evolu- 
tionary tracks for pre-main sequence stars 
from Siess et al. (2000). We also over- 
plot the position of Cernis 52 which seems 
to be consistent with a theoretical 2.0 ± 
0.2 M Q star with an age of 3-20 Myr. The 
position of the star is in between the theo- 
retical tracks with masses 2.2 and 1.8 M . 
These masses and the derived log g imply 
a radius in the range 1.1 < R/R Q < 2.8. 

The young cluster IC 348 has an age 
of 3-7 Myr, according to Trullols & Jordi 
(1997), from a sample of 123 stars in a wide 
range of spectral types. Their analysis re- 
lies on photometric data and is possibly 
affected by unknown extinctions. Luhman 
et al. (2003) discuss the age of IC 348 in 
detail and claimed that the observational 



data seems to imply ages ranging from 1 
to 10 Myr, although it favours a mean age 
of 2 Myr. They used photometric data of a 
large sample of K7-M8 spectral-type stars, 
although they also obtained spectroscopic 
data. Herbig (1998) also studied the age of 
roughly 100 stars of this cluster and found 
a spread between 0.7 and 12 Myr. All 
these studies depend on spectral-type clas- 
sifications which may be uncertain, how- 
ever, they seem to support that the age of 
IC 348 likely range between 1 and 12 Myr. 
The position of Cernis 52 in Fig. [10] is in 
agreement with the most likely age range 
for the cluster IC 348. 

7.3. Luminosity and Distance 

There is no available paralax for Cer- 
nis 52. Here, we estimate the distance 
to Cernis 52 from different photomet- 
ric magnitudes taking into account the 
stellar parameters. We know the photo- 
metric magnitudes in five different filters: 
my = 11.4±0.1 (from Cernis 1993), m# = 
10.65 ± 0.10 and mj = 10.04 ± 0.10 mag, 
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and mj = 8.91 ± 0.02, m H = 8.49 ± 0.03, 
and rriK B = 8.27 ± 0.02 mag from the 
2MASS catalog^. Except for the 2MASS 
magnitudes, we adopted ad hoc the values 
for the uncertainties on the photometric 
magnitudes of Cernis 52. We derived the 
radius of the star from the surface gravity, 
log(g/cm s 2 ) = 4.2±0.4 dex, and assuming 
a mass of 2 M . This radius, together with 
the spectroscopic estimate of the effective 
temperature, T e s = 8350 ±200 K, provides 
an intrinsic bolometric luminosity in the 
range 5.4 < L+/L Q < 34.2. We determined 
the absolute bolometric magnitude from 
the following formula: Mbd,* — Mbd,© = 
—2.5 log(L*/ Lq), where the solar bolomet- 
ric magnitude is = 4.75 (IAU sim- 
posium 1999) and the solar luminosity is 
L Q = 3.847 x 10 33 erg s -1 . 

The apparent bolometric magnitude, 
m t>oi,* = m i,o ~ BCi, was derived using 
the bolometric corrections, BCi, for non- 
overshooting ATLAS 9 models (Bessell, 
Castelli, & Plez 1998). The bolometric 
corrections were determined as BCi = 
BC v + V-i. For the 2MASS infrared mag- 
nitudes we computed the theoretical colors 
V — i, for i = J, H, K s , from Gonzalez 
Hernandez & Bonifacio (2009). 

We adopt the color excess, E(B — V) = 
0.9 mag, from Cernis (1993) which seems 
to provide consistent distances when us- 
ing optical and infrared magnitudes. We 
compute the magnitude corrected for ex- 
tinction in each filter as m^o = rrii — Ai, 
where Ai were obtained using the relation 
Ai = RiE(B - V), with R iy the ratio of to- 
tal to selective extinction, given by the co- 
efficients provided in McCall (2004). This 



3 http : / / www. ipac . caltech.edu / 2mass/ 



author gives Ry = 3.07 with deviations un- 
likely to exceed 0.05. Note that using a 
different value, e.g. Ry = 3.31 as derived 
for Cernis 77 (BD+31° 643) in Snow et al. 
(1994) or R v = 3.2 as in Cernis (1993), 
lead to small corrections to the distance 
determination, by -21 and -12 pc respec- 
tively. 

The apparent magnitude, rrii, was de- 
contaminated from the magnitude of the 
companion star, according to the stellar 
parameters derived in Section |H We esti- 
mate the flux ratio of both stars, F sec / Fc52, 
from theoretical magnitudes which were in- 
terpolated in the grid of theoretical colors 

V — R and V — I (Bessel et al. 1998) and 
the color V — K and magnitudes J, H, K s 
(Gonzalez Hernandez & Bonifacio 2009). 
Finally, the corrected magnitude is derived 
from the following expression 

mi = m 2 i,i + log(F 2M + 1) 

where rrii is the corrected magnitude in 
the filter i, rri2i,i is the observed magnitud 
which contains the flux of both stars, and 
i^i,!, the flux ratio F sec /Fc52 in the filter 
%. In Table |3] we show the derived distance 
for each filter with different sets of the rel- 
evant parameters. Note that the relatively 
small formal error in the distance estimate 
for each filter is calculated by assuming the 
magnitudes equal rrii + Am;. 

We have derived a distance of 23llg 35 pc 
for the star Cernis 52, according to the 
adopted stellar parameters and mass in 
section [5\2l and I7T21 Cernis (1993) argued 
that the star Cernis 52 is a member of the 
young cluster IC 348. This author derived 
a distance to Cernis 52 of 236 pc from the 

V magnitude, adopting a spectral type of 
A3V and a color excess E(B — V) — 0.9. 
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Fig. 10. — Theoretical evolutionary tracks 
from Siess et al. (2000) for pre-main se- 
quence stars with masses 1.8, 2.0 and 2.2 
M and for solar metallicity, Z = 0.02. 
The different ages of the tracks refer to the 
age of a star at the end of the track, i.e. 
close the main sequence. The filled circle 
represents the position of Cernis 52. The 
triangles and squares following the tracks 
represent the position of the theoretical 
stars at the age of 3 and 9 Myr respec- 
tively. 

He also derived a distance to the cluster 
IC 348 of 260 ± 16 pc from 13 probable 
members of the cluster. Trullols & Jordi 
(1997) provide a photometric determina- 
tion of the distance to IC 348 based on 
the filter V of 240lgf pc from 43 clus- 
ter members. These distance determina- 
tions rely on spectral type classification. 
Herbig (1998) argued that spectral types 
derived from spectroscopy would produce 
types somewhat later, and thus smaller 
Ay, and concluded that the Trullols & 
Jordi (1997) photometric distance should 
correspond to a spectroscopic distance of 
276 pc. Our distance determination to 
Cernis 52 seems to be in agreement with 



that of the young cluster IC 348. 

To sum up, all the arguments provided 
in this paper support Cernis 52 as being a 
member of the young cluster IC 348, and 
place the star far enough to explain the 
large color excess, and the presence of the 
interstellar features recently discovered by 
Iglesias-Groth et al. (2008). These intestel- 
lar features must be caused by the known 
dark clouds already mentioned by Cernis 
(1993), in particular, the dark cloud L1470 
which covers all the young cluster IC 348. 

7.4. Photospheric abundances 

The O I triplet 6156-8 A required a high 
O abundance ([0/H]= 0.5), whereas the 
best fit abundance of the O I triplet 7771- 
5 A is [O/H] = 0.15 (see Fig.E}. This is un- 
expected since the O I triplet at 7771-5 A 
shows stronger negative NLTE corrections 
than the O I triplet 6156-8 A (Baschek 
et al. 1977). This may indicate the exis- 
tence of some veiling at 7773 A, instead of 
the zero value adopted (according to the 
extrapolation of the veiling previously de- 
rived function of wavelength) . Similar 
problems are found in the Mg lb triplet 
5167-83 A (see Fig. Ej). For these features 
we depict the synthetic spectra computed 
using the best fit abundance given by the 
feature Mg I 5173 A, and this abundance 
is unable to reproduce both Mg features at 
5167 and 5183 A. There are other features 
over the whole spectrum that cannot be 
reproduced with the synthetic stellar spec- 
tra which may be related to the presence 
of unknown DIBs. We emphasize that for 
the chemical abundance analysis we have 
chosen only those features that are not af- 
fected by any known interstellar absorp- 
tion band. 
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The photospheric abundances of the 
star seems to be almost solar within the 
error bars, i.e., the star does not belong 
to the group of chemically peculiar (CP) 
A-type stars. There appears to be a cor- 
relation between the presence of chemi- 
cal anomalies and the rotational velocity 
of these stars (e.g. Fossati et al. 2008; 
Takeda et al. 2008). Thus, stars with 
v sin i < 50 — 75 km s _1 are underabundant 

in O and Ca (down to [X/H] 0.8), and 

overabundant in Fe ([Fe/H]~ 0.5) and spe- 
cially Ba (up to [Ba/H]~ 1.5), whereas Si 
and S shows near-solar abundances. Cer- 
nis 52 has a rotational velocity on the edge 
between "peculiar" and "normal" A stars. 
In the bottom panel of Fig. [TTJ we de- 
pict the Ba line which gives an abundance 
[Ba/H]~ —0.8, which is even lower than 
the solar abundance. This definitively dis- 
cards Cernis 52 SIS db chemically peculiar 
A-type star. 

7.5. The naphthalene feature at 
6707.4 A 

In our paper announcing the discovery 
of the naphthalene cation (Iglesias-Groth 
et al. 2008), we noted that the strongest 
band of the interstellar cation appeared 
at 6707.4 A, i.e. very close to the stel- 
lar Li I resonance line. The abundance 
analysis presented in this paper was ex- 
pected to provide indications of the likely 
Li abundance for the star, i.e., a star with- 
out abundance anomalies is very likely to 
have a lithium abundance no greater than 
the local abundance of A(Li)= 3.3 (see 
bottom panel of Fig. [9]). 

In Fig. [12] we show three synthetic spec- 
tra: (i) two of them computed only consid- 
ering the star Cernis 52 and with a already 




6096 6098 6]00 6102 6104 6106 6108 6110 
A. (ANGSTROMS) 




6135 6140 6145 6150 

A. (ANGSTROMS! 



Fig. 11. — Same as Fig. [TJ but for other 
spectral regions. Upper panel: synthetic 
spectra are computed with A(Li) = 4.0 
(solid line) and A(Li) = 3.3 (dashed line). 
We emphasize that A(Li) = 3.3 must be 
the maximum Li abundance in this star 
(see Section [73]) . Bottom panel: synthetic 
spectra are computed with [Ba/H] = —0.8 
(solid line) and [Ba/H] = —0.1 (dashed 
line) . 

too high Li abundance, A(Li) = 4.0 and 
a cosmic Li abundance, A(Li) = 3.3, and, 
(ii) the other combining the spectrum of 
Cernis 52 with A(Li) = 3.3 and the spec- 
trum of the companion star with A (Li) = 
4.0, contributing with 90 and 10%, respec- 
tively, to the stellar flux. We estimated 
the rotational velocity (see Section [3]) at 
vsini ~ 65 km s _1 from many stellar lines 
in the spectral range AA5270 — 6400 A. As 
seen in Fig. [9] and more clearly in Fig. [121 
the stellar lithium line at 6707.8 A is un- 
able to fill the whole feature even if we 
increase the Li abundance up to A (Li) = 
4.0. In the second case, we choose a ro- 
tational velocity for the companion star 
of 100 km s _1 to reproduce this broad fea- 
ture, but even using a very high Li abun- 
dance, A (Li) = 4.0, is not enough to com- 
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Fig. 12. — Same as Fig. but for other 
spectral region. Synthetic spectra are com- 
puted with A (Li) = 4.0 (solid line) and 
A(Li) = 3.3 (dashed-dotted line). We also 
display the combined spectrum (dashed 
line) including the spectrum of Cernis 52 
contributing with 90% of the stellar flux 
and the spectrum of the companion star 
contributing at 10%, respectively. The 
dashed line shows a synthetic spectrum 
with A(Li) = 3.3 and 4.0, for Cernis 52 and 
the companion, respectively. 

pletely fill this broad feature. The re- 
quired Li abundance able to fit this feature, 
A(Li) = 5.6, is extraordinarily high. Again 
A(Li) = 3.3 is a highly probably maximum 
for the companion at T eff = 5800 K. We 
must emphasize that even with this large 
rotational velocity and abnormally high Li 
abundance for the companion star, the fea- 
ture is still not perfectly reproduced by the 
synthetic spectrum. 

In the upper panel of Fig. [LTJ, we de- 
pict the subordinate Li I line at 6103.4 A. 
The subordinate Li line is very weak and 
the feature at 6103 A is dominated by a 
Ca I line at 6102.2 A and a Fe II line at 
6103.5 A. Thus, increasing the Li abun- 
dance from A(Li)= 3.3 to A(Li) = 4.0 



has a negligible effect on the line profile. 
Thus, we cannot estimate the Li abun- 
dance from this subordinate Li line. Polo- 
sukhina & Shavrina (2007) suggested that 
the Li abundance obtained in CP stars 
with magnetic fields from the Li subordi- 
nate line may be higher than that derived 
from the resonance Li line by 0.2-0.4 dex. 
This may suggest vertical stratification of 
Li. In addition, CP stars shows enhanced 
Li abudances up to A(Li) = 4.0. 

In Section 17.41 we discarded the possi- 
bility that Cernis 52 is a CP star, so its Li 
abundance must be at the most the cosmic 
Li abundance, i.e. A(Li) = 3.3. Therefore, 
as seen in Fig. [TJl the stellar Li line should 
only be responsible for a small fraction of 
the equivalent width of the broad feature 
at 6707.4 A. 

All the above statements confirm that 
this broad feature must have a insterstel- 
lar origin and most probably linked to the 
naphthalene molecule. The other broad 
feature at 6756 A (see Fig. [9]) is proba- 
bly another interstellar feature (S. Iglesias- 
Groth et al. 2009, in preparation). 

Finally, another detection of PAHs in 
the interstellar medium has been recently 
discovered by Iglesias-Groth et al. (2009). 
These authors report a PAH band at ~ 
7088 A, the anthracene cation, CuKf . 

8. Conclusions 

We have performed a detailed chemical 
analysis of the star Cernis 52. We apply 
a technique that provides a determination 
of the stellar parameters, taking into ac- 
count any possible source of veiling. We 
find T eff = 8350 ± 200 K, log(#/cm s 2 ) = 
4.2±0.4, [Fe/H] = -0.01±0.15, andaveil- 
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ing (defined as F vei i/F tota i) of less than 55% 
at 5000 A and decreasing toward longer 
wavelengths. 

The spectrum of Cernis 52 shows many 
features very likely related with the inter- 
stellar medium. In addition, we discover in 
photometric images a companion 1.7 mag 
fainter star at a distance of 0.818"±0.007", 
but this star only contributes with 10% of 
the stellar flux in our spectra and hardly 
affect the stellar features of Cernis 52. The 
derived chemical abundances are roughly 
solar within their error bars. This prevent 
the star from being a chemically peculiar 
star. 

We have determined the radial velocity 
of Cernis 52 at v r = +13.7 ± 1 kins" 1 , 
being almost equal to the mean radial ve- 
locity of the young cluster IC 348. 

We have compared the stellar parame- 
ters with pre-main-sequence evolutionary 
tracks of solar metallicity and see that 
the star is consistent with being a pre- 
main-sequence A-type star with an age of 
3-20 Myr. 

We have estimated the distance to 
Cernis 52 using the available Johnson- 
Cousins and 2MASS photometric data, 
at 23llgg 5 pc according to the stellar pa- 
rameters and its error bars. This value 
also agrees with the distance to the cluster 
IC 348. 

The proper motion of Cernis 52, (/i a cos S, 
fi s = (+7.22,-8.62) ± (1.5,1.6) mas yr" 1 , 
is consistent with the proper motion of 
IC 348. 

All these measurements make it likely 
that the star Cernis 52 (BD+31° 640) be- 
longs to the young cluster IC 348. 

We confirm that the feature at 6707.4 A 



is not related with a stellar lithium line 
because the rotational velocity of the 
star, vsini = 65 + 5 kms -1 , is too low 
to explain the broad feature associated 
with the naphthalene cation. Further- 
more, the presence of a companion star 
cannot either explain this feature even 
for an abnormally high Li abundance, 
log[iV(Li)/iV(H)] + 12 > 5 dex. 

As already stated in Cernis (1993), 
the interstellar features, in particular, the 
naphthalene cation, that appear in the 
spectrum of Cernis 52 may form in the 
dark cloud L1470 which covers all the clus- 
ter IC 348 and is at about the same dis- 
tance. 

J.I.G.H. acknowledges support from 
the EU contract MEXT-CT-2004-014265 
(CIFIST). D.L.L. thanks the Robert A. 
Welch Foundation of Houston, Texas, for 
support through grant F-634. D.A.G.H. 
acknowledges support from the Span- 
ish Ministry of Science and Innovation 
(MICINN) under the 2008 Juan de la 
Cierva Programme. This work has been 
partially founded by project AYA2007- 
64748 of the Spanish Ministry of Education 
and Science. We thank the FastCam team 
at Instituto de Astrofisica de Canarias and 
Universidad Politecnica de Cartagena for 
kindly obtaining images of Cernis 52 for 
this work. We are grateful to Victor J. 
S. Bejar and Jose A. Caballero for help- 
ful discussions. We are grateful to Tom 
Marsh for the use of the MOLLY analy- 
sis package. This work has made use of 
the VALD, SIMBAD, DSS and 2MASS 
databases, and IRAF and Aladin facilities. 
The Two Micron All Sky Survey is a joint 
project of the University of Massachusetts 



19 



and the Infrared Processing and Analy- 
sis Center /California Institute of Technol- 
ogy, funded by the National Aeronautics 
and Space Administration (NASA) and 
the National Science Foundation (NSF). 
The compressed files of the "Palomar Ob- 
servatory - Space Telescope Science Insti- 
tute Digital Sky Survey" of the northern 
sky, based on scans of the Second Palomar 
Sky Survey are copyright (c) 1993-2000 by 
the California Institute of Technology and 
are distributed herein by agreement. All 
Rights Reserved. Produced under Con- 
tract No. NAS5-2555 with the National 
Aeronautics and Space Administration. 

REFERENCES 

Al-Naimiy, H. M. 1978, Ap&SS, 420, 183 

Ali, A. W., & Griem, H. R. 1965, Phys. 
Review, 140, 1044 

Ali, A. W., & Griem, H. R. 1966, Phys. 
Review, 144, 366 

Andrievsky, S. M., Chernyshova, I. V., 
Usenko, I. A., Kovtyukh, V. V., 
Panchuk, V. E., & Galazutdinov, G. A. 
1995, PASP, 107, 219 

Barklem, P. S., Stempels, H. C, Allende 
Prieto, C, Kochukhov, O. P., Piskunov, 
N., & O'Mara, B. J. 2002, A&A, 385, 
951 

Baschek, B., Scholz, M., & Sedlmayr, E. 
1977, A&A, 55, 375 

Bessell, M. S., Castelli, F., & Plez, B. 1998, 
A&A, 333, 231B 

Castelli, F., & Kurucz, R. L. 2001, A&A, 
372, 260 



Cernis, K. 1993, Baltic Astronomy, 2, 214 

Cowley, C. R., & Castelli, F. 2002, A&A, 
387, 595 

Dahm, S. E. 2008, AJ, 136, 521 

Draine, B. T., & Lazarian, A. 1998, ApJ, 
494, L19 

Ecuvillon, A., Israelian, G., & Santos, N. 
C. 2006, A&A, 445, 633 

Evans, D. S. 1967, Determination of Radial 
Velocities and their Applications, 30, 57 

Fossati, L., Bagnulo, S., Landstreet, J., 
Wade, G., Kochukhov, O., Monier, R., 
Weiss, W., & Gebran, M. 2008, A&A, 
483, 891 

Gonzalez Hernandez, J. I., et al. 2004, 
ApJ, 609, 988 

Gonzalez Hernandez, J. I., Rebolo, R., 
Israelian, G., Casares, J., Maeda, K., 
Bonifacio, P., & Molaro, P. 2005, ApJ, 
630, 495 

Gonzalez Hernandez, J. I., Rebolo, R., & 
Israelian, G. 2008a, A&A, 478, 203 

Gonzalez Hernandez, J. I., Rebolo, R., 
Israelian, G., Filippenko, A. V., 
Chornock, R., Tominaga, N., Umeda, 
H., & Nomoto, K. 2008b, ApJ, 679, 732 

Gonzalez Hernandez, J. I., & Bonifacio, P. 
2009, A&A, 497, 497 

Gratton, R. G., et al. 2001, Experimental 
Astronomy, 12, 107 

Grevesse, N., Noels, A., & Sauval, A. J. 
1996, The sixth annual October Astro- 
physics Conference, ASP Conf. Ser. 99, 
117 



20 



Griem, H. R. 1960, ApJ, 132, 883 

Herbig, G. H. 1998, ApJ, 497, 736 

Iglesias-Groth, S., Manchado, A., Garcia- 
Hernandez, D. A., Gonzalez Hernandez, 
J. L, & Lambert, D. L. 2008, ApJ, 685, 
L55 

Iglesias-Groth, S., Manchado, A., Rebolo, 
R., Gonzalez Hernandez, J. I., Garcfa- 
Hernandez, D. A., & Lambert, D. L. 
2009, ApJ, submitted 

Kurucz, R. L. ATLAS9 Stellar Atmo- 
spheres Programs and 2 km s _1 Grid 
(CD-ROM, Smithsonian Astrophysical 
Observatory, Cambridge, 1993) 

Kurucz, R. L. 2005, Memorie della Societa 
Astronomica Italiana Supplement, 8, 14 

Kurucz, R. L., Furenild, I., Brault, J., & 
Testerman, L. 1984, Solar Flux Atlas 
from 296 to 1300 nm, NOAO Atlas 1 
(Cambridge: Harvard Univ. Press) 

Li, A. 2009, Deep Impact World Ob- 
servatory Event: Synergies in Space, 
Time, and Wavelength, 161 

Loktin, A. V., & Beshenov, G. V. 2003, 
Astronomy Reports, 47, 6 

Luhman, K. L., Stauffer, J. R., Muench, 
A. A., Rieke, G. H., Lada, E. A., Bou- 
vier, J., & Lada, C. J. 2003, ApJ, 593, 
1093 

McCall, M. L. 2004, AJ, 128, 2144 

Nordhagen, S., Herbst, W., Rhode, K. L., 
& Williams, E. C. 2006, AJ, 132, 1555 

Oscoz, A., et al. 2008, Proc. SPIE, 7014 



Piskunov, N. E., Kupka, F., Ryabchikova, 
T. A., Weiss, W. W., & Jeffery, C. S. 
1995, A&AS, 112, 525 

Polosukhina, N. S., & Shavrina, A. V. 
2007, Astrophysics, 50, 381 

Roser, S., Schilbach, E., Schwan, H., 
Kharchenko, N. V., Piskunov, A. E., & 
Scholz, R.-D. 2008, A&A, 488, 401 

Rosolowsky, E. W., Pineda, J. E., Foster, 
J. B., Borkin, M. A., Kauffmann, J., 
Casein, P., Myers, P. C, & Goodman, 
A. A. 2008, ApJS, 175, 509 

Sbordone, L. 2005, Memorie della Societa 
Astronomica Italiana Supplement, 8, 61 

Scholz, R.-D., et al. 1999, A&AS, 137, 305 

Siess, L., Dufour, E., & Forestini, M. 2000, 
A&A, 358, 593 

Sneden, C. 1973, PhD Dissertation, Univ. 
of Texas, Austin 

Snow, T. P., Hanson, M. M., Seab, G. C, 
& Saken, J. M. 1994, ApJ, 420, 632 

Takeda, Y., Han, I., Kang, D.-L, Lee, B.- 
C, & Kim, K.-M. 2008, Journal of Ko- 
rean Astronomical Society, 41, 83 

Trullols, E., & Jordi, C. 1997, A&A, 324, 
549 

Tull, R. G., MacQueen, P. J., Sneden, C, 
Lambert, D. L. 1995, PASP, 107, 251 

Watson, R. A., Rebolo, R., Rubiho- 
Martm, J. A., et al. 2005, ApJ, 624, L89 

Watson, R., Rebolo, R., Davis, R. J., & 
Rubino-artm, J. A. 2006, CMB and 
Physics of the Early Universe, 20-22 
April 2006, Ischia, Italy, p.66 



21 



This 2-column preprint was prepared with the 
AAS IATfiX macros v5.2. 



22 



